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Abstract 

Q-i We investigate the prospects for detecting and measuring the parameters of 



WIMP dark matter in a model independent way at the International Linear Collider. 
The signal under study is direct WIMP pair production with associated initial state 
radiation e + e~ —> XXI- The analysis accounts for the beam energy spectrum of the 
ILC and the dominant machine induced backgrounds. The influence of the detector 
parameters are incorporated by full simulation and event reconstruction within the 
framework of the ILD detector concept. We show that by using polarised beams, 
the detection potential is significantly increased by reduction of the dominant SM 
background of radiative neutrino production e + e~ — > vPj. The dominant sources 
of systematic uncertainty are the precision of the polarisation measurement and the 
shape of the beam energy spectrum. With an integrated luminosity of C = 500 fb _1 
the helicity structure of the interaction involved can be inferred, and the masses 
and cross-sections can be measured with a relative accuracy of the order of 1%. 



To be submitted to EPJC 



1 Introduction 



Weakly interacting massive particles (WIMPs) \ with masses in the order of M x ~ 
100 GeV are among the favoured candidate particles to provide the observed cosmo logical 
abundance of dark matter. Complementary to direct and indirect detection experiments 
which look for signals of primordial WIMPs, colliders could produce WIMP particles 
under laboratory conditions. Since the WIMPs themselves leave collider experiments 
undetected due to their weak interaction with matter, collider searches typically rely on 
WIMPs appearing in cascade decays of more heavy exotic particles, thus assuming a 
specific extension of the Standard Model of particle physics (SM). 

Direct WIMP production does not depend on the existence or kinematical accessibility of 
additional new particles. However it can be observed only when accompanied by wide- 
angle initial state radiation (ISR). It has been shown pQ that for annihilation cross-sections 
compatible with the relic abundance of dark matter, the cross-section for radiative WIMP 
production at an electron-positron Linear Collider can be sufficiently large to observe this 
process above the irreducible background from radiative neutrino production. The energy 
spectrum of the ISR photon can be exploited to extract information on the WIMP mass 
and cross-section. The resulting discovery reach and mass resolutions have been studied 
before assuming either supersymmetric extension of the SM or the presence of universal 
extra-dimensions [2]. In both cases the new partners of the electron can be exchanged 
in the t-channel and thus impact the shape of the photon energy spectrum. Based on 
four-vector smearing and taking into account the radiative neutrino background and an 
overall systematic uncertainty of 0.3%, the two cases can be distinguished for a large 
range of the parameter space, and the mass of the electron partner can be measured far 
beyond the kinematic limit for its direct production. 

However, the shape of the photon spectra from signal and SM background are expected to 
be significantly distorted by various effects related to the detector, the reconstruction and 
the accelerator. The observation reach for the single photon signature at the International 
Linear Collider (ILC) has been investigated in full detector simulation in [3J. Considering 
the statistical uncertainty only and assuming fully polarised beams, it has been shown 
that with an integrated luminosity of 500 fb" 1 cross-sections down to about 12 fb can be 
observed at the 5cr level. 

This publication studies the precision to which the WIMP mass and cross-section can 
be determined based on full detector simulation and taking into account all relevant 
backgrounds as well as systematic uncertainties of the detector measurement and the 
limited knowledge of beam parameters. Instead of studying different explicit extensions 
of the SM, the generic parametrisation of the WIMP cross-section pQ is used here to 
investigate if the dominant partial wave of the WIMP production cross-section can be 
determined. 

Beyond the WIMP approach, the results of this study are relevant for pair production of 
any invisible or nearly invisible particle, whenever the ISR recoil method is applicable. 

This paper is organised as follows: In the next section, the model-indepedent Ansatz for 
radiative WIMP production as well as the signal charcteristics and the most important 
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backgrounds are introduced. Section [3] discusses the experimental conditions at the ILC 
and one of the proposed detector concepts, including the resulting event selection criteria 
and systematic uncertainties. The achievable precisions on the WIMP mass and its po- 
larised and unpolarised cross-sections as well as the prospects to determine the dominant 
partial wave of the WIMP production mechanism are finally presented in Section |H 



2 Radiative WIMP production in e + e collisions 



WIMPs can be pair produced in e + e~ collisions if they have a non-vanishing coupling to 
electrons and if their mass doesn't exceed half of the center-of-mass energy. An additional 
ISR photon allows to detect such events. In this analysis no specific scenario of physics 
beyond the SM is assumed and thus only radiation off the incoming particles is considered, 
as illustrated in Figure [TJ 

Figure 1: Illustration of the 
radiative WIMP pair production 
mechanism. No explicit assumption 
on the WIMP-fermion interaction 
is made, thus only ISR off the in- 
coming particles is considered. 



The production cross-section for WIMP pairs with an associated ISR photon of energy 
Ej and polar angle 9 can be written in the limit of non-relativistic final state WIMPs 
as PQ: 

da CWCeffaal + fl-z)* 2Jn , n „ ( . N 1/2+J ° 




dx d cos 6 167T xsin 2 # 



Here, M x is the WIMP mass, S x its spin, s the center-of-mass energy squared and the 
dimensionless variable x = Jo is the quantum number of the dominant partial wave 
in the production process. The quantitiy K e is the annihilation fraction of WIMPs into 
electron positron paira^. It implicitely depends on the helicity structure of the WIMP 
production mechanism and the helicities of the beam electrons and positrons. Under 
the additional assumption that the primordial dark matter consists of our WIMPs, the 
overall scale of the production cross-section above is given by the WIMP annihilation 
cross-section into fermion-antifermion pairs a an . It is estimated from the observed relic 
density to be about 0.8 pb for s-wave and about 7 pb for p-wave annihilation [I]. 

Figure [2] shows the total unpolarised cross-section for e + e~ — > XXI a t a center-of-mass 
energy of 500 GeV as a function of the WIMP mass for p-wave production and a (a) 
spin-1/2 or (b) spin-1 WIMP. The phase space of the photon has been restricted to an 

1 With the sum over all final state fermions (i = e, u, q, ...) J2i K i = 1- 
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Figure 2: Total visible cross-sections for radiative WIMP production at 500 GeV for 
K e = 1: (a) p-wave, spin- 1/2, (b) p-wave, spin-1. The error bars illustrate an 10% 
uncertainty. 



experimentally accessible range of > 8 GeV and |cos#| < 0.995. The error bars 
illustrate a 10 % uncertainty, the annihilation cross-sections have been taken from pQ. 
The cross-section is well above 100 fb for WIMP masses up to 180 GeV in the spin- 1/2 
and up to 200 GeV in the spin-1 case. For K e in the order of 10% and large WIMP masses, 
the expected cross-sections get close to the minimum of 12 fb for a 5 a observation with 
an integrated luminosity of 500 fb~ x found in a previous study [3J. 

While the angular distribution of the ISR photons is independent of the following hard 
process, the photon energy spectrum strongly depends on the WIMP parameters. Fig- 
ure |3ja) shows the shape of the photon energy distributions expected for different WIMP 
masses in the case of p-wave production, featuring the mass dependent cut-off at the max- 
imal allowed photon energy. The shape of the photon energy threshold near its endpoint, 
which is equivalent to the rise of the production cross-section at threshold, is compared in 
Figure [3Jb) for the cases of J = (s-wave) and Jo — 1 (p-wave) production. Measurement 
of the threshold behaviour would provide an indication of the partial wave involved in the 
production process. This information could for instance be employed to test whether the 
WIMP is a majorana fermion (like for instance a neutralino) or a scalar particle. 

The single photon signature suffers from substantial irreducible background from radiative 
neutrino production e + e~ — > vvj, which has an unpolarised cross-section in the order of 
several pb, depending on the acceptance cuts on the photon energy and angle. Due to 
the size of this background, a similar analysis at LEP would not have been feasible. At 
the ILC, the much higher luminosity and the beam polarisation help to gain sensitivity 
over the neutrino background, which proceeds primarily via t-channel W exchange and 
hence is strongly polarisation dependent. Other important background processes comprise 
multi-photon final states e + e~ — > 77 and radiative Bhabha scattering e + e~ — > e + e~7, as 
well as machine induced backgrounds. 
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Figure 3: Generator level photon energy spectra a) for different WIMP masses in the 
p-wave case b) for M x = 160 GeV in p- and s-wave production. All spectra are normalised 
to an integral of 1 to facilitate the shape comparison. 
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Process 


(-0.8; +0.3) 


(+0.0; +0.0) 


(+0.8; -0.3) 




5821 


2575 


1263 




782.0 


355.4 


214 




55.8 


26.2 


19 


77 




11.4 x 10 3 




777 




1.1 x 10 3 




7777 




0.1 x 10 3 




e + e~ 




890 x 10 3 





Table 1: Cross-sections of the most important Standard Model backgrounds for three 
different polarisation configurations. All photons denoted explicitely in the process name 
are included in the matrix element. In addition two ISR photons are included. At least 
one photon is required to have P 7 > 8 GeV and | cos6*| < 0.995. 



Table [T] shows the cross-sections of these Standard Model backgrounds for three different 
beam polarisation configurations as obtained from Whizard |4]. P e -and P e + denote the 
values of the electron and positron beam polarisation, respectively. All photons denoted 
explicitely in the process name are included in the matrix element calculation performed 
by O'Mega [5]. At least one photon is required to be in the analysis phase space as defined 
above. In each case two additional soft ISR photons are included in the event generation. 

The comparison of these numbers to the expected signal cross-sections displayed in Fig- 
ure [2] immediately demonstrates the need for a careful evaluation of the experimental 
possibilities to suppress all reducible background processes. Furthermore the systematic 
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uncertainties associated with the remaining irreducible background need to be understood 
in order to obtain a realistic picture of the ILC's capabilities in this channel. 

Since on an event-by-event basis the radiative WIMP production is indistinguishable from 
radiative neutrino production, the signal sample was generated by reweighting one third 
of the vi?'-f sample according to the ratio of the tree-level differential cross-sections of 
WIMP and neutrino production. This method has the benefit of giving flexible access 
to the all combinations of WIMP properties without the need for dedicated Monte-Carlo 
production at each considered point in parameter space. 



3 Experimental Conditions 



This study is based on an integrated luminosity of 500 ftr 1 and the nominal parameter 
set of the ILC at a center-of-mass energy of 500 GeV as specified in its Reference Design 
Report [5J. For comparison, the SB2009 beam parameter set [7] is used for selected 
results. The chosen beam parameter sets differ mainly in the strength of the focussing of 
the beam at the interaction point, resulting in different beam energy spectra, which have 
been derived using GuineaPig [8]. 

The electron beam polarisation of \P e -\ = 80% and the positron polarisation of \P e +\ = 
30% are assumed to be known to 5P/P = 0.25% In order to estimate the impact 
of these quantities, alternative values of \P e +\ = 60% and 5P/P = 0.1% are used for 
comparison. The peak beam energy will be measured to a relative precision of 10~ 4 [H] 
by the energy spectrometers. For the evaluation of the systematic uncertainties in this 
analysis, a slightly more conservative estimate based on the position of the radiative 
return to the Z in the photon energy spectrum is used. Assuming the same amount of 
integrated luminosity as used in the main analysis, this position can be determined to a 
statistical precision of 80 MeV. For the integrated luminosity a relative precision of 10 -3 
is assumed [TO] . 

The analysis is performed in full simulation of the ILD detector concept [11] . The central 
device for detecting the ISR photons is the electromagnetic calormeter (ECAL), which is 
a highly granular SiW sampling calorimeter with a cell size of 5 x 5 mm 2 . In testbeam 
measurements [12j, a resolution of 



M = (16 . 6±0 .i)% x ^e(i.i ± o.i)%, (2) 

E \JE 

as has been achieved. The simulation used here shows a similar performance. 

ECAL clusters to which no track can be associated are considered as photon candidates. 
After corrections for split clusters and losses in the cracks between modules, the energy 
of isolated photons can be reconstructed without bias over the whole polar angle range. 
This is illustrated by Figure Ufa), which shows the ratio of reconstructed over true photon 
energy as function of cos# of the photon. Residual deviations from unity near module 
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Figure 4- (a) average of reconstructed over true photon energy as function of cos 9 after 
corrections for cluster splitting and losses in insensitive regions of the detector, 
(b) reconstruction efficiency for isolated photons as function of cos 9. 



boundaries are smaller than 2%. Figure ID^b) shows the efficiency of the photon recon- 
struction as a function of cos 9. The average efficiency is near 97%, varying from close to 
98% in the middle of the detector to about 95% in the endcaps. For further analysis, only 
events with at least one photon with a reconstructed energy 10 GeV < E™ c < 220 GeV 
and a polar angle | cos# rcc | < 0.98 are considered. 

In order to distinguish photons from electrons or hadrons and to veto any other significant 
detector activity, the tracking system and the scintillator-tile steel sandwich hadronic 
calorimeter (AHCAL) are employed in addition. Both ECAL and AHCAL as well as the 
tracking system are placed inside a superconducting coil which provides a 3.5 T solenoidal 
magnetic field. The tracking system combines a Time Projection Chamber (TPC) with 
additional silicon pixel and strip detectors to achieve a tracking efficiency of > 99.5% in 
an angular range of 7° < 9 < 173° (| cos#| < 0.9925). The veto on tracks of charged 
particles needs to be adjusted to allow for possible tracks from beam background or low 
p t photon-photon interactions. 

The transverse momentum distribution of these tracks obtained from full detector simula- 
tion and reconstruction is displayed in Figure O together with the corresponding distribu- 
tion for tracks originating from e + e~ pair background. Both distributions are normalised 
to the number of tracks expected from these processes per bunch crossing (BX), i.e. 0.7 
tracks/BX from 77 processes and 1.5 tracks/BX from e + e~ pairs. In the analysis, events 
containing any track with p? > 3 GeV are vetoed. The remaining probability to falsely 
reject a signal event due to an overlayed background track is about 0.25%. In addition to 
the track veto, events with other neutral particles are rejected by requiring that the total 
visible energy of the event must not exceed the recontructed photon energy by more than 
20 GeV. 
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Figure 5: Transverse momen- 
tum distribution of tracks from e + e~ 
pair background and from 77 inter- 
actions. The distributions are nor- 
malised to the number of tracks ex- 
pected per bunch crossing from these 
sources. 
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The low angle calorimeters complement the hermeticity of the detector, leaving only 
holes of 4.5 mrad and 5.6 mrad in the polar angle coverage of the BeamCal around 
the incoming and outgoing beams, respectively. The innermost part of the BeamCal is 
exposed to energy depositions of several TeV per bunch crossing from e + e~ pairs created 
by Beamstrahlung processes. Isolated high energy electrons (or photons) can be detected 
above this background with efficiencies larger than 80% for angles above typically 20 mrad 
from the detector z axig^] [13]. The complete angular and energy dependence of this 
efficiency has been taken into account in the simulation. This is crucial in order to obtain 
a correct estimate of the background from radiative Bhabha scattering, since vetoing 
such isolated high energy electrons in the BeamCal is mandatory to reduce the Bhabha 
background to a manageable level. 



Process 


E^; cos 9 


PT,track 


E ■ - E 

-'-'vis J - /r f 


BeamCal tag 


Eff. 


%] 






2493.3 


2435.4 


2283.88 


2239.63 


8 


9.5 




i/z/77 


344.3 


325.4 


238.52 


228.51 


66.4 


Z/P777 


25.4 


23.2 


11.82 


11.05 


43.5 


77 


578.1 


457.3 


60.74 


5.80 




1.0 


777 


145.0 


112.7 


4.65 


0.10 




0.1 


7777 


19.5 


14.7 


0.15 


0.03 




0.2 


e + e~ 


421533.1 


88935.9 


67389.80 


1228.70 




0.3 



Table 2: Number of events for the main background processes at the different selection 
stages for an integrated luminosity of C — 1 fb" 1 and (P e -;P e +) = (+0.0; +0.0). The 
second column contains the event numbers in the analysis phace space. The third column 
lists the event numbers after the cuts on the energy and polar angle of the reconstructed 
photon. The last column lists the selection efficiencies. 

2 The detector z axis points in the middle between incoming and outgoing beams, which cross under 
an angle of 14 mrad. 
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Table |2] illustrates the effect of the selection cuts for the case of unpolarised beams. Beam 
polarisation significantly reduces or enhances the rates for the neutrino processes while 
leaving the Bhabha rate at the same level. The final selection efficiency for the missing 
energy plus photon signature is near 90%. In real data, this efficiency could be controlled 
from the high energy part of the photon energy spectrum which contains the radiative 
return to the Z and is not used in the WIMP analysis as well as from radiative muon pair 
production u + u~j. From the statistics available in these two processes, we estimate the 
systematic uncertainty on the selection efficiency to 1.5%. In addition, the impact of the 
cut on the reconstructed photon energy depends slightly on the WIMP mass. Based on 
the accuracies on the WIMP mass achieved in this study, this translates into an additional 
sytematic effect on the selection efficiency. Added in quadrature, both effects together 
yield a total systematic uncertainty of 1.75% on the selection efficiency. 




V] 



Figure 6: Photon energy spectra after all selection cuts for Standard Model background 
and a 150 GeV p-wave WIMP signal, normalised to an integrated luminosity of 50 fb^ 1 and 
beam polarisations of (P e -; P e +) = (+0.8; —0.3). (a) Background (red), signal (green) and 
a statistically independent "data" sample (points) (b) Signal and "data" after background 
subtraction. 

The photon energy spectrum obtained after all selection cuts is displayed in Figure [B£a) 
for Standard Model background (red) and for a WIMP signal (green; M x = 150 GeV, 
p-wave). The distributions are normalised to an integrated luminosity of 50 fb _1 and 
beam polarisations of (P e -;P e +) = (+0.8; —0.3). The "data" points are obtained from a 
statistically independent Monte-Carlo sample assuming the presence of background and 
signal. Figure M)o) shows the "data" after subtracting the expected background compared 
to the signal expectation (green). The size of the statistical errors of the "data" points 
reflects the fluctuations of the large amount of subtracted background. 
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4 Results 



In this section the precisions on the polarised and unpolarised cross-sections and on the 
WIMP mass achievable at the ILC are presented. Additionally the possibility to dis- 
criminate between s- and p-wave production is investigated. In the model-independent 
approach, the branching fraction of WIMP annihilation into e + e~ pairs, K e , has a free 
dependency on the helicity of the initial state electrons. In order to illustrate the power 
of polarised beams to determine this dependency, three different coupling scenarios are 
compared in the following: 

• "Equal": The WIMP couplings are independent of the helicity of the incoming 
electrons and positrons, i.e. «;(e^,e^) = «;(e^,e^) = K(e^,e£) = /c(e^,e^). 

• " Helicity" : The couplings conserve helicity and parity, 

K ( e Ri 4) = «(ez, 4); 4) = K ( e Z> 4) = °- 

• "Anti-SM": This scenario is a "best case" scenario, since the WIMPs couple 
only to right-handed electrons and left-handed positrons: «(e^,ej), with all other 
n(e-,e + ) = 0. 

4.1 Cross-section Measurements 

For this first part of the analysis a typical running scenario of the ILC is assumed, where an 
integrated luminosity of C = 500 fb _1 is split into 200 fb" 1 each with (+|P e - 1; — \P e +\) and 
(— |P e - |; + |P e +|) (short (H — ) and ( — h)), as well as into 50 fb -1 each with (+|P e -|; +|P e +|) 
and (-|P e -|; -|P e +|) (short (++) and ( )). 

By combining the four independent measurements with different beam polarisation con- 
figurations, the helicity structure of the WIMP interactions can be extracted by solving 
the equation system 

P e -\)(l + \P e+ \)a LL 
P e -\)(l-\P e+ \)a LR } 
P e -\)(l-\P e+ \)a LL 
P e -\){l + \P e+ \)a LR } 
P e -\)(l-\P e+ \)a LL 
P e -\)(l + \P e+ \)a LR ] 
P e -\)(l + \P e+ \)a LL 
P e -\)(l-\P e+ \)a LR }, (3) 

In practice, \P e -\ (as well as |P e +|) will not be exactly equal for the four running periods, 
so that the measured cross-sections need to be extrapolated to the same |P e -| (|P e +|) 
based on the actual polarimeter measurements and their systematic uncertainty. 
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Figure 7: Fully polarised cross-sections o"{j?,l} measured within three WIMP scenarios: 
(a) + (b) "Equal", (c) + (d) "Helicity", (e) + (f)' "Anti-SM", and for two different absolute 
values of the positron polarisation: (a) + (c) + (e) \P e +\ = 30%, (b) + (d) + (f) \P e +\ = 60%. 
The cross-sections are normalised to the input unpolarised cross-section ctq. 
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(|P e -|;|P e+ |) = (a 


.8; 0.3) 


(|P e -|;|P e+ |) = (0 
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0.99 ± 0.24 


(0.16) 
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(0.15) 


"Helicity" scenario 


vrl/vo 


1.99 ± 0.24 
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(0.08) 


vrr/vq 
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3.99 ±0.12 


(0.10) 


vrr/vq 


0.00 ±0.33 


(0.22) 


0.00 ±0.23 


(0.14) 


vll/(?o 
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(0.28) 


0.00 ±0.23 


(0.15) 


0~LR./o~o 
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Table 3: Fully polarised cross-sections o"{_r,l} measured within three WIMP scenarios 
and for two different absolute polarisations of electrons and positrons. The values are 
normalised to the input unpolarised cross-section of ao = 100 fb. The quoted uncertainties 
are the squared sum of statistical errors and systematic uncertainties, with the bracketed 
values corresponding to an increased precision on the polarisation measurement ofSP/P = 
0.1%. 

The results for the three studied coupling scenarios and a WIMP mass of 150 GeV are 
listed in terms of the fully polarised cross-sections G{r,l\ i n Tableland shown in Fig- 
ures[7](a) to (f). Independently of the coupling scenario, the precision on the fully polarised 
cross-sections is dominated by the uncertainty on the beam polarisation as it enters the 
extrapolation to \P e -\ = \P e - \ = 100%. This can be seen from the fact that the uncer- 
tainty reduces considerably when the positron polarisation is increased from 30% to 60% 
(right column in Table |3J). About the same improvement is obtained by reducing the 
uncertainty of the polarisation measurement from 0.25% to 0.1% (values in parentheses 
in Table ED- In any of the cases, the coupling scenarios can be clearly separated from each 
other. 

The unpolarised cross-section cro can be measured by combining all four measurements 
without the need to extrapolate to \P e -\ = \P e -\ = 100%, which reduces impact of the 
polarimeter uncertainties significantly. Table H] gives the achievable precisions on o~q in 
the three different coupling scenarios. Even though its impact is reduced, the polarisation 
uncertainty still dominates the total error of 4 to 5 fb for SP/P = 0.25% , which can be 
considerably reduced to 2 to 3 fb when assuming 5P/P = 0.1%. 
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Figure 8: Sensitivity in terms of standard deviations a as a function of the unpolarised 
WIMP cross-section for the three coupling scenarios and different assumptions on the 
positron polarisation and on the precision of the polarisation measurement. 

13 



Data scenario Unpolarised cross-section: cr 


-1- ctat" 
Zl olcLL 


4- 


sys 


(± total) [fb] 


(simulated) ( | P e - | ; P e + ) 


= (0.8; 0.3) 


CI P 


~ 1 1 


ln + |) 


= (0.; 


3; 0.6) 


Assumed polarisation uncertainty 


oPj P = O.zo/o 














11 TTi ill (\ f\ r\ i o o 1 

Equal 99.0 ± 2.8 ± 


4.3 (± 5.1) 


QQ 9 


4- 


2.7 


± 


3.5 


(± 4.4) 


Helicity 99.1 ± 2.3 ± 


4.0 (± 4.6) 


99.4 


± 


2.0 


± 


2.8 


(± 3.4) 


Anti-olVl 99. o ± 1.4 ± 


2.o (± 3.2) 


99.7 


± 


1.1 


± 


O 1 


/_|_ o A\ 

(± 2.4) 


Assumed polarisation uncertainty 


= 0.10% 














"Equal" 99.0 ± 2.6 ± 


2.0 (± 3.3) 


99.1 


± 


2.6 


± 


1.9 


(± 3.2) 


"Helicity" 99.1 ± 2.3 ± 


2.0 (± 3.0) 


99.3 


± 


2.0 


± 


1.8 


(± 2.6) 


"Anti-SM" 99.6 ± 1.4 ± 


1.8 (±2.3) 


99.7 


± 


1.2 


± 


1.7 


(± 2.1) 



Table 4 : Measured unpolarised cross-section by a combination of cross-section mea- 
surements with polarised beams for an integrated luminosity of C = 500 fb -1 . 



These numbers have been derived based on an unpolarised cross-section of cro = 100 fb. 
However the measurements are dominated by systematic uncertainties and by the statis- 
tical uncertainty on the background. Therefore the derived uncertainties are to a large 
extent independent of the actual signal cross-section. Figure [8] shows the sensitivity in 
terms of standard deviations (a) as a function of the unpolarised WIMP cross-section 
for the three coupling scenarios and different assumptions on the positron polarisation 
and on the precision of the polarisation measurement. In absence of a signal in the data, 
cross-sections above 8.2 (3.6) fb can be excluded at 90% CL in the worst (best) case. 
Cross-sections above 25 (11) fb could be observed at the 5a level in the worst (best) case, 
which is in reasonable agreement with the minimal observable cross-section found in [3] 
based on statistical uncertainties only. 

4.2 Mass Measurement and Extraction of the Dominant Partial 
Wave 

For the measurement of the WIMP mass and the extraction of the dominant partial wave, 
an integrated luminosity of 500 fb" 1 of the (|P e -|; — \P e + 1) configuration is assume ci. The 
WIMP mass is determined by comparing template photon energy spectra for background 
and signals of different masses against the "data" spectrum and by searching the mass for 
which the x 2 is minimised. This method is much more powerful than a mere determination 
of the maximal photon energy, because near the endpoint of the spectrum the signal is 
buried in the fluctuations of the huge Standard Model background, as can be seen in 
Figure [6]^b). 

3 This is effectively an additional 300 fb -1 with respect to the data set assumed for the cross-section 
measurement, which corresponds to a little more than one year of operation time of the ILC after the 
initial four years. 
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Figure 9: Reconstructed vs true WIMP mass for the right and wrong partial wave 
assumptions in the templates, (a) true s-wave, (b) true p-wave. 



The shape of the photon energy spectrum below the endpoint however also depends on 
the dominant partial wave of the production, c.f. Figure [3](b) . Figure [9] shows the impact 
of the partial wave assumption in the templates on the fitted WIMP mass. While the 
reconstructed mass follows the true mass nicely when the correct partial wave templates 
are used, the wrong templates yield fitted masses either 20 GeV too low ((a), true s-wave) 
or too high ((b), true p-wave). 

However the correct assumption on the dominant partial wave can be easily determined 
over the full WIMP mass range by comparing the x 2 P er degree of freedom of the template 
fits if positron polarisation is available. This can be seen in Figure [TQl where the x 2 / n df 
is displayed as a function of the true WIMP mass for the right and wrong partial wave 
hypothesis. In the upper row, the positron polarisation is set to \P e +\ = 0, and in this 
case no convincing difference in the \ 2 values is obtained. For \P e +\ = 30% (middle row) 
and \P e + \ = 60% (bottom row), the wrong assumption yields a x 2 / n df significantly larger 
than unity. Since for this test it is crucial to have a simulated sample corresponding to 
the actually assumed intergrated luminosity, the limited amount of events available from 
full simulation has been augmented by a fast simulation sample obtained from an ILD 
version of SGV [H] for this (and only this) figure. 

After resolving the ambiguity due to the dominant partial wave, the main systematic 
uncertainty in the mass measurement is the shape of the beam energy spectrum followed 
by the uncertainty on the beam energy scale. The former has been estimated from a 
comparison to fits with templates based on the beam energy spectrum obtained with 
the SB2009 beam parameter set. Conservatively, the full resulting shift in measured 
mass is included in the systematic errors quoted here. Thus the uncertainty of the mass 
measurement can be reduced by determining the beam parameters, for instance from the 
shape of the energy depositions of the e + e~ pair background in the forward detectors [13] . 
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Figure 10: y 2 /ndf of the template fit using the right and wrong partial wave assumptions 
as a function of the true WIMP mass. (a)+(b): \P e +\ = 0, (c)+(d): \P e +\ = 30%, (e) + (f): 

\P e +\= 60% 



16 



(P /P_) = (0.8/-0.3) 



J = 1 ; Equal 




-0.05 



100 120 140 160 180 200 220 240 

IVL [GeV] 



(P /P D ) = (0.8/-0.3) 



J = 1 ; Helicity 




-0.05 



100 120 140 160 180 200 220 240 



IVL [GeV] 



;0.05 
>0.04 
3 0.03 
0.02 
0.01 


-0.01 
-0.02 
-0.03 
-0.04 



-0.05 



(P/Pp) = (0.8/-0.3) 
J = 1 ; Anti-SM 




n ... I ... I ... I ... I ... I ... I ... I 



100 120 140 160 180 200 220 240 



IVL [GeV] 



(P e /P p ) = (0.8/-0.6) 
J = 1 ; Equal 




-0.05 



100 120 140 160 180 200 220 240 

M x [GeV] 



(P e /Pp) = (0.8/-0.6) 
J = 1 ; Helicity 




-0.03 
-0.04 
-0.05 



(d) 



100 120 140 160 180 200 220 240 

M x [GeV] 



;0.05 
>0.04 
3 0.03 
0.02 
0.01 


-0.01 
-0.02 
-0.03 
-0.04 
-0.05 



(P e /P p ) = (0.8/-0.6) 

J = 1 ; Anti-SM 




■i ... i ... i ... i ... i ... i ... i ... i 



(f) 



100 120 140 160 180 200 220 240 

M x [GeV] 



Figure 1 1 : Relative uncertainty on the reconstructed WIMP mass as a function of the 
true WIMP mass for the three different coupling scenarios and two different values of the 
positron polarisation. The blue area shows the systematic uncertainty and the red bands 
the additional statistical contribution. 
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Mass 


WIMP mass: ±stat. 


± SE (scale) ± 5E (shape) (total) [GeV] 


[GeV] 


(P e -; P e+ ) = (0.8; -0.3) 


(P e -; P e+ ) = (0.8; 


-0.6) 


" Equal" scenario 


120 


2.48 ±0.07 ±1.90 (3.12) 


2.24 ±0.07 ±1.90 


(2.93) 


150 


1.98 ±0.05 ± 1.46 (2.46) 


1.83 ±0.05 ±1.45 


(2.33) 


180 


1.69 ±0.03 ± 1.00 (1.96) 


1.57 ± 0.03 ± 1.00 


(1.86) 


210 


0.79 ±0.02 ±0.54 (0.96) 


0.74 ±0.02 ±0.54 


(0.91) 


"Helicity" 


scenario 






120 


1.92 ± 0.07 ± 1.89 (2.70) 


1.53 ±0.07 ±1.89 


(2.43) 


150 


1.62 ±0.05 ±1.46 (2.18) 


1.23 ±0.05 ±1.45 


(1.90) 


180 


1.36 ±0.03 ± 1.00 (1.69) 


0.94 ± 0.03 ± 1.00 


(1.37) 


210 


0.67 ±0.02 ±0.54 (0.87) 


0.59 ±0.02 ±0.54 


(0.80) 


"Anti-SM" 


scenario 






120 


1.04 ± 0.07 ± 1.88 (2.15) 


0.82 ±0.07 ±1.88 


(2.05) 


150 


0.81 ±0.05 ± 1.45 (1.66) 


0.72 ±0.05 ±1.44 


(1.61) 


180 


0.66 ±0.03 ±1.00 (1.19) 


0.37 ± 0.03 ± 1.00 


(1.06) 


210 


0.16 ±0.02 ±0.55 (0.59) 


0.09 ±0.02 ±0.55 


(0.59) 



Table 5: Statistical and systematic uncertainties on the measured WIMP masses for an 
integrated luminosity of C = 500 fb _ in the three coupling scenarios "Equal", "Helic- 
ity" and "Anti-SM" for three different polarisation configurations. 



This is forseen in the design of the ILC detectors, but has not been studied in this 
analysis. Table and Figure [TT] display the statistical and systematic accuracies reached 
under these assumptions for the three coupling scenarios and two values of the positron 
polarisation. The total uncertainty ranges between 0.5 and 3%. For higher WIMP masses 
the measurement becomes in general more precise, because the photon energy spectrum 
is restricted to lower values and thus the signal is more prominent locally. 

5 Conclusions 

The photon plus missing energy signature at the International Linear Collider has been 
investigated in the context of a model-independent characterisation of WIMPs. In view of 
the substantial Standard Model backgrounds to this signature, the analysis has been per- 
formed in full detector simulation and with realistic assumptions on the beam properties 
and resulting systematic uncertainties. The unpolarised cross-section can be measured 
with an accuracy between 2 and 5 fb depending on the beam polarisation and the WIMP 
scenario. The helicity structure of the WIMP-fermion interaction can be obtained from 
the polarised cross-section. These measurements are systematically limited by the knowl- 
edge of the beam polarisation. 
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Via the shape of the photon energy spectrum, p- and s-wave production can be clearly 
distinguished if positron polarisation is available. Finally the mass of the WIMP can be 
determined with an accuracy between 0.5% and 3%, depending on the WIMP parameters. 
These numbers are limited by the knowledge of the shape of the beam energy spectrum. 

Beyond the model-independent WIMP approach, these results can also be used in any 
specific scenario beyond the Standard Model with invisible or nearly invisible massive 
particles in the kinematic reach of the ILC. 
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